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We present the first study of the time-dependent CP- violating asymmetry in B — ► J/ip tt decays 
using e + e~ annihilation data collected with the BABAR detector at the T(4S) resonance during 
the years f 999-2002 at the PEP-II asymmetric-energy B Factory at SLAC. Using approximately 88 
million BB pairs, our results for the coefficients of the cosine and sine terms of the CP asymmetry 
are C J/i}7V o = 0.38 ± 0.41 (stat) ± 0.09 (syst) and S J/i}7T o = 0.05 ± 0.49 (stat) ± 0.16 (syst). 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 

The Standard Model of electroweak interactions de- plex phase in the three-generation Cabibbo-Kobayashi- 
scribes CP violation in B-mcson decays by a com- Maskawa (CKM) quark-mixing matrix [1]. The b — > ccs 
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modes such as B° — > J/tpK® yield precise measurements 
of the quantity sin2,3, where (3 = &rg{—V cd V* b /V td V t * b ] 
(see for example Refs. [2-4]). The decay B° — ► J/tpn° is 
a Cabibbo-suppressed b — > ccd transition. In the Stan- 
dard Model both £?° -» J/tpK° and 5° -» J/V-tt have 
penguin amplitudes with the same weak phase as the 
tree amplitude, and an additional penguin amplitude 
with a different phase. In B° — ► J/ipK®, the penguin 
amplitude with a different weak phase is suppressed by 
X CKM , where Xckm is the sine of the Cabibbo angle, 
while in B° — > J/tpn , the tree and each penguin am- 
plitude are equal to leading order in Xckm- Therefore, 
B° — ► J/tpn may have a CP asymmetry that differs 
from that of B° — ► J/tpK®, with the size of the asymme- 
try serving as a probe of the penguin decay amplitudes 
in both modes. 

BABAR has previously measured the B° — ► J/tpn° 
branching fraction L (2.0±0.6 (stat)±0.2 (syst))xl(T 5 [5], 
using T(4S) — ► SB decays. For the CP asymmetry mea- 
surement, the flavor (£? or B°) of the B meson that 
decays to J/tpn is inferred, or tagged, using properties 
of the other B meson and the time evolution of the BB 
system. The decay time distributions, f + (f_), of B de- 
cays to a CP cigenstate with a B° (B°) flavor tag, are 
given by 



f±(Af) 



-\At\/r B o 

4r s o 



liSj/^o sin(Am d At) 



TCj/^o cos (Am d At) 



(1) 



where At = i roc — t tag is the difference between the proper 
decay time of the reconstructed B meson and the proper 
decay time of the tagging B meson, t b o is the B° life- 
time, and Anid is the B°-B° oscillation frequency. The 
coefficients can be expressed in terms of a complex pa- 
rameter A, which depends on both the B°-B° oscilla- 
tion amplitude and the B° and B° decay amplitudes 
to this final state [6]: Sj/^^o = 2TmA/(l + |A| 2 ) and 
Cj/Vtt" = (1 - |A| 2 )/(1 + |A| 2 ). A decay amplitude with 
only a tree component would give Sj/^^a — — sin2/3 and 



C 



= 0. 



The data used in this measurement were collected with 
the BABAR detector [7] at the PEP-II storage ring in the 
years 1999 to 2002. Approximately 81 fb _1 of e + e~ anni- 
hilation data recorded at the T(4S) resonance are used, 
corresponding to a sample of approximately 88 million 
BB pairs. An additional 5fb _1 of data collected approx- 
imately 40MeV below the T(45) resonance are used to 
characterize non-BB background sources. 

B° — > J/tp 7T° candidates are selected (details are given 
in Ref. [5]) by identifying J/tp — ^ e + e~ or J/tp — ► \x + \x~ 



decays and 7r° — > 77 decays. For the J/tp 



TABLE I: Efficiencies for the requirement on the Fisher dis- 
criminant and flavor tagging, given independently, with sta- 
tistical uncertainties. 



Type of event 


Efficiency (%) 




Fisher 


Tagging 


B° -> J/V>7T° 


99.2 ± 0.1 


65.6 ± 0.6 


B° -> J/iP K° s (ir Tr°) bkg. 


98.9 ± 0.1 


65.6 ± 0.6 


Inclusive J/tp bkg. 


94.9 ± 0.7 


70.4 ± 1.4 


BB generic bkg. 


98.5 ± 0.4 


61.1 ± 1.6 


Continuum bkg. 


28.6 ± 0.7 


52.3 ± 0.8 



mass of the lepton pair is required to be between 2.95 and 
3.14GeV/c 2 , and 3.06 and 3.14GeV/c 2 , for the electron 
and muon channels, respectively. The photon candidates 
used to reconstruct the n° candidate are identified as 
clusters in the electromagnetic calorimeter (EMC) with 
polar angles between 0.410 and 2.409 rad, that are spa- 
tially separated from every charged track, and have a 
minimum energy of 30 MeV. The lateral energy distribu- 
tion in the cluster is required to be consistent with that 
of a photon. The invariant mass of the photon pair is re- 
quired to between 100 and 160MeV/c 2 . Finally, the J/tp 
and 7T° candidates are assigned their nominal masses and 
combined using 4-momentum addition. 

Two kinematic consistency requirements are applied 
to each B candidate. The difference, AE, between the 
£?-candidate energy and the beam energy in the e + e~ 
center-of-mass (CM) frame must be —0.4 < AE < 
0.4 GeV. The beam-energy-substitutcd mass, mss = 
must be greater than 5.2 GeV/ c 2 , 



— > fi + /1 ) channel, each lepton candidate must be consis- 
tent with the electron (muon) hypothesis. The invariant 



where y/s is the total CM energy and p* B is the B- 
candidate momentum in the CM frame. 

A linear combination of several kinematic and topo- 
logical variables, determined with a Fisher discriminant, 
provides additional separation between signal and e + e~ 
—> QQ (<Z = u , d, s, c) continuum background events. The 
Fisher discriminant uses the following inputs: the zeroth- 
and second-order Legendre polynomial momentum mo- 
ments (L = |p*| and L 2 = £ 4 \p*\ *£°£°i=± , where 
p* are the CM momenta for the tracks and neutral 
calorimeter clusters that are not associated with the sig- 
nal candidate, and 6^ are the angles between p* and 
the thrust axis of the signal candidate); the ratio of 
the second-order to zeroth-order Fox- Wolfram moments, 
again using just tracks and clusters not associated with 
the signal candidate; | cos 0t \ , where Qt is the angle be- 
tween the thrust axis of the B candidate and the thrust 
axis of the remaining tracks and clusters in the event; 
and I cos Oe | , where 0g is defined as the angle between the 
negative lepton and B candidate directions in the J/tp 
rest frame. The requirement placed on the Fisher dis- 
criminant is 99% efficient for signal and rejects 71% of 
the continuum background. The efficiencies for satisfy- 
ing this requirement are summarized in Table I. 
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We split the backgrounds into four mutually exclu- 
sive categories, two of which have a J/ip from B de- 
cays (B — ► J/ip X). The first background category is 
B° — > J/ip K^(tt tt ) decays where one of the ir° mesons 
is nearly at rest in the e + e~ CM frame. The second 
background category consists of other B — > J/ip X de- 
cays (inclusive J/ip), which contribute through random 
combinations of J/ip and 7T° candidates. The third and 
fourth categories consist of random combinations of par- 
ticles in BB decays (BB generic) and continuum events, 
respectively. Monte Carlo simulation [8] is used to model 
aspects of the B° — > J/ip K®(ir°ir ), inclusive J/ip, and 
BB generic backgrounds. A sample (J/V'fakc) selected 
from data taken below the T(4S) resonance is used to 
model the continuum background. In this case, the J/ip 
candidate is reconstructed from two tracks that are not 
consistent with a lepton hypothesis. Monte Carlo simula- 
tion is used to check that this procedure, which increases 
the size of the sample, correctly models the continuum 
background. 

The algorithm for £>-flavor tagging assigns events to 
one of four hierarchical, mutually exclusive tagging cat- 
egories, and is described in detail in Ref. [3]. The to- 
tal tagging efficiency for the signal and each background 
source is given in Table I. Untagged events are excluded 
from further consideration. Vertex reconstruction and 
the determination of At follow the techniques detailed in 
Ref. [9]. We require —20 < At < 20 ps and an estimated 
uncertainty on At of less than 2.4 ps. 

We extract the CP asymmetry by performing an un- 
binned extended maximum likelihood fit. The likeli- 
hood is constructed from the probability density func- 
tions (PDFs) for the variables toes, AE, and At. The 
quantity that is maximized is the logarithm of 



C = 



N 5 

IIE 

i=l j=l 



(2) 



where 



is the number of events for each of the five 



hypotheses (one signal and four background) and N is 
the number of input events. The Vj are the one- or two- 
dimensional PDFs for variables d, for each signal or back- 
ground type. The parameters /"* are the tagging frac- 
tions for each of the tagging categories on (assigned for 
each event i) and each of the signal or background types 
j. For the B° -> J/iPir° signal and B° -> J/ip K°(ir°ir°) 
background, the values of /"* are measured with a sam- 
ple (-Bflav) of neutral B decays to flavor eigenstates con- 
sisting of the channels D^*^h + (h + = n + ,p + , and af) 
and J/iPK*°(K*° -> K+tt~) [3]. Monte Carlo simulation 
is used to estimate the values for the inclusive J/ip 
and BB generic backgrounds, while the •^/V'fake 

sample 

is used for the continuum background. 

The signal toes distribution is modeled as the sum 
of two components. The first is a modified Gaussian 
function that, for values less than the mean, has a 



width parameter that scales linearly with the distance 
from the mean. The second component, accounting for 
less than 6% of the distribution, is a threshold func- 
tion [10], which is a phase-space distribution of the form 

toes^/(1 ~ |P^) exp(£(l - fjp^)), with a kinematic 
cut-off at -Ebcam = 5.289 GeV and one free parameter £. 
The signal AE distribution is modeled by the sum of 
a Gaussian core with an asymmetric power-low tail [11] 
and a second order polynomial. The parameters of these 
PDFs are determined by fitting to a signal Monte Carlo 
sample. The peak position of the AE distribution is a 
free parameter of the full CP likelihood fit to allow for 
EMC energy scale uncertainties. 

The kinematic variables toes and AE are corre- 
lated in the B° -» J/ip K° s (tt tt°) and inclusive J/ip 
backgrounds, so two-dimensional PDFs are employed 
for these modes. Variably-binned interpolated two- 
dimensional histograms of these variables are constructed 
from the relevant Monte Carlo samples. 

The toes PDFs for the BB generic and continuum 
backgrounds arc modeled by the threshold function given 
above, and the AE PDFs for these two backgrounds are 
modeled by second order polynomials. The parameters 
for these PDFs are obtained from the BB generic Monte 
Carlo sample and the J/ip {ake sample. 

The PDFs used to describe the At distributions of 
the signal and background sources are each a convolu- 
tion of a resolution function 1Z and decay time distribu- 
tion V: V(At,a At ) = K(6t,a A t) ® £>(Ai true ), where At 
and Attruo are the measured and true decay time differ- 
ences, St = At — Attruo, and o At is the estimated event- 
by-event error on At. 

For the signal, the resolution function consists of the 
sum of three Gaussian distributions, the parameters of 
which are determined from the -Bflav sample, as in the 
B° — ► J/ip K® measurement [9] . The decay time distribu- 
tion is given by Eq. 1 modified for the effects of _B-flavor 
tagging: 

p -|At|/r B o 

V± J (At)=—^-— {(l T A Wa ) 

± S f (1 - 2w a ) sin(ATO d At) 
T Cf (1 - 2w a ) cos(Am d At)}, 



(3) 



where T^Xf^af) i s f° r a B°(B°) tagging meson. The 
variable w a is the average probability of incorrectly tag- 
ging a B° as a B° or a B° as a B° (wf°), and 
Aw a — — . Both w a and Aw„ are deter- 
mined using the Ba av data sample [3] . We use the values 
Am d = 0.489 ps" 1 and t b o = 1.542 ps [12]. 

The PDF used to model the At distribution for the 
B° — > J/ipK^(TT°Tr°) background, which also includes a 
CP asymmetry, takes the same form as that for signal, 



but with S 



j/4,k° = sin2/3 = 0.74 [3] and C J/4 , K o 



0. 



The parameterizations of the At PDFs for the inclu- 
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TABLE II: Results of the CP likelihood fit, for the full region 
-0.4 < AE < 0.4 GeV and m ES > 5.2GeV/c 2 . Errors are 
statistical only. The global correlation coefficient is 0.14 for 
Cj/^o and 0.15 for S J/4 , n o. 



Fit results 




0.38 


± 


0.41 


Sj/i>n° 


0.05 


± 


0.49 


Signal AE peak position ( MeV) 


-13.2 


± 


7.2 


B° — > J/ip-K signal (events) 


40 


± 


7 


B° -> J/ipKK^-K ) background (events) 


140 


± 


19 


Inclusive J/tj) background (events) 


109 


± 


35 


BB generic background (events) 


52 


± 


25 


Continuum background (events) 


97 


± 


22 



sive J/ip and BB generic backgrounds each consist of 
prompt and exponential decay components. Decays ap- 
pear to be prompt when particles from the reconstructed 
B are erroneously included in the tagging B vertex. For 
the BB generic background, the prompt and exponen- 
tial components correspond to the cases where the two 
decay products forming the J/ip come from both or just 
one of the B mesons, respectively. The fraction that is 
in the exponential component, the decay lifetime param- 
eter, and the resolution parameters are determined from 
the Monte Carlo simulation. 

The At PDF for the continuum background has only a 
prompt component and the resolution parameter values 
are obtained by fitting the J/ - 0f a ke sam pl e - 

The results of the CP asymmetry fit, for all free pa- 
rameters, are shown in Table II. There are 40 ± 7 signal 
events in the total sample of 438 selected events. The 
projection in ibes is shown in Fig. 1. The yields and 
asymmetry as functions of Ai, overlaid with projections 
of the likelihood fit results, are shown in Fig. 2. Repeat- 
ing the fit with the added constraint Cj/^^o = does not 
significantly change the result for Sj/^ T o . 

The dominant contributions to the systematic errors 
in Cj/^^o and Sj/^^o are summarized in Table III. The 
first class of uncertainties are those obtained by varia- 
tion of the parameters used in the /ties, AE, and At 
PDFs, where the dominant sources are the uncertain- 
ties in the signal AE PDF parameters. A systematic 
error to account for a correlation between the tails of 
the signal thes and AE distributions is obtained by us- 
ing a two-dimensional PDF. Another contribution stems 
from the impact of EMC energy scale uncertainties on 
the modeling of the B° — > J/tp K^(Tr°Tr°) background. 
An additional systematic uncertainty comes from the 
choice of the binning of the two-dimensional PDFs for 
the B° — > J/tp Kg(TT°Tr ) and inclusive J/tp backgrounds. 

In summary, an unbinncd extended maximum like- 
lihood fit yields 40 ± 7 signal events and the param- 
eters of time-dependent CP asymmetry for the decay 
B° -> J/iPn°: Cj/^o = 0.38 ± 0.41 (stat) ± 0.09 (syst) 
and Sj/^o = 0.05 ±0.49 (stat) ±0.16 (syst). Within the 
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FIG. 1: Projection in mES for the results of the CP fit, dis- 
played with the added requirement —0.11 < AE < 0.11 GeV. 
In contrast, the CP fit uses the full AE region. In the further 
restricted region tties > 5.27 GeV/ c 2 , there are 49 data events 
(points), of which about 12 events are fit as background. Here, 
B° — ► J/ip Kg(ir°n ) and inclusive J/ip decays contribute to 
the enhancement in the background distribution at large mES • 
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FIG. 2: Distributions of events a) with a B° tag (N B o), b) 
with a B° tag (Ng ), and c) the raw asymmetry (N B o — 
Ngo)/(N B o + Ngo), as functions of At. Candidates in these 
plots are required to satisfy —0.11 < AE < 0.11 GeV and 
m-ES > 5.27 GeV/ c 2 . Of the 49 signal and background events 
in this region, 25 have a B° tag and 24 have a B° tag, with 
fit background contributions of approximately 5 and 7 events, 
respectively. The curves are projections that use the values 
of the other variables in the likelihood to determine the con- 
tributions to the signal and backgrounds. 



Standard Model formulation of CP asymmetries, these 
results demonstrate the possibility, with additional inte- 
grated luminosity, of observing penguin contributions in 
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TABLE III: Summary of systematic uncertainties. 



Source 




Sj/ij, n o 


Parameter variations 


m-ES and AE parameters 


0.05 


0.13 


Tagging fractions 


0.00 


0.01 


At parameters 


0.03 


0.02 


Additional systematics 


A_E-m.ES correlation in signal 


0.07 


0.08 


EMC energy scale B° -> J/ip K%(n°iT ) 


0.01 


0.00 


Choice of two-D histogram PDFs 


0.01 


0.03 


Beam spot, boost /vtx., misalignment 


0.01 


0.01 


Total systematic uncertainty 


0.09 


0.16 



B a — > J/tpir°. Such a measurement may experimentally 
constrain similar amplitudes in B° — ► J/ipK°. 
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